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Abstract: We evaluate a new 3-loop sum-integral which contributes to the Debye screening 
mass in hot QCD. While we manage to derive all divergences analytically, its finite part is 
mapped onto simple integrals and evaluated numerically. 
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1 Introduction 

Motivated by the need to compute a number of unknown three-loop vacuum sum-integrals of 
mass dimension 2, in order to complete the evaluation of the 3-loop Debye screening mass 
in hot QCD [1], we turn our attention towards one particular open case. The present paper 
is hence rather technical in practice, but serves to establish one of the few missing building 
blocks for the determination of the Debye mass at NNLO. 

In Ref. [1], the bosonic contribution to has been reduced to a sum of ten sum-integrals 
times pre-factors which are rational functions in d, by systematic [2] use of integration- 
by-parts (IBP) relations [3], generalized to finite temperatures [4]. Two out of the basis 
of those ten master sum-integrals were triple products of 1-loop cases I and hence trivial 
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(cf. Eq. (B.2)), while only two further non-trivial three-loop representatives of that basis 
(containing basketball-type sum-integrals) are presently known, namely S\ of [5] (see also 
[6, 7]) and -B32 of [1], however not deeply enough in their e-expansion since the corresponding 
pre- factors develop single poles as d— >3. 

It can be shown that, by a change of basis, the eight non-trivial three-loop basketball- 
type masters can be traded for only three spectacles-type (plus some trivial factorized) sum- 
integrals, with pole-free pre-factors [8]. Of those three new masters, it turns out that one has 
already appeared in a study of the g 7 contributions to the pressure of massless thermal scalar 
</> 4 theory in the framework of screened perturbation theory. It has originally been computed 
in Ref. [9] and been re-derived in Ref. [10], where it was named Mifi. This leaves us with 
two elements of the new basis to be evaluated, one of which we tackle in the present paper. 

Let us define the massless bosonic 3-loop vacuum sum-integral V2 in terms of the 1-loop 
2-point sum-integrals II, II as 

We use (Euclidean) bosonic four-momenta P = (Pq, p) = (2-Kn p T, p) with P 2 = P$ + p 2 ; the 
temperature of the system is denoted by T; and the sum-integral symbol is a shorthand for 

^ r £/(Sp. witW = 3-2e. (1.2) 

In the remainder of this paper, we evaluate the new sum-integral V2 analytically as far as 
possible, utilizing methods from [5, 7, 10, 11]. In Sec. 2 we split up the sum-integral into 
nine pieces that are either calculable analytically, or explicitly finite such that they can be 
evaluated numerically in d = 3. Sec. 3 treats each of these nine terms in turn, which are then 
summed to obtain our final result in Sec. 4. We conclude in Sec. 5, while some technicalities 
are relegated to the appendices. 

2 Decomposition of V 2 

Inspired by [11] and guided by our experience from [10], V2 can be identically re-written as 

v, = f p n(P) fi(P) + $-^{ P - n B ] [n - n B ] + 
+ u D [n - Be] + [n B - n D ] [n - n c ] + [n - n c ] u D + [n - n c ] [fl B - u D ] + 
+ [n c - n B ] n B + n B [n c - n B ] + n B n B } , (2.1) 

= Vi + V 2 + V 3 + V 4 + V 5 + V 6 + V 7 + V 8 + V 9 , (2.2) 

where 5p picks out the Matsubara zero-mode, the primed sum excludes the zero-mode and 
we have suppressed the argument (P) of all functions in curly brackets. Our strategy amounts 
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to choosing (IIb: zero-T part of II; [lie — II5]: large-P behavior of remainder) 



1 



G(l,l,d + 1) 



Q Q 2 (P-Q) 2 (p2 )( 3-d)/2 > 



n c = n B + 



27? 

p2 ' 



nD - Pf^ ' (2 - 3) 



where the functions G and I are known analytically as given in App. B while a\ (and a% 
below) are constants to be fixed later, as well as (with U = (1,0)) 



^0 -u u 

R R 2 {P-R) 2 ~ At v 



RnR U 



A(P 2 ) + PiB{P 2 ) = { — + 1)P ° 2 ^ 



R R 2 (P-R) 2 



U„U V {g^A(P 2 ) + P, t P u B(P 2 )} 



p2 to , 9 r2 



Ad 

{d+l)P 2 -P 2 G(l,l,d + 1) 



P 2 4d (a 2 r 2 )( 3 - rf )/ 2 

3 Evaluation of the nine terms of Eq. (2.1) 



(2.4) 
(2.5) 



We will now turn to each term of Eq. (2.1) in sequence, evaluating it up to the constant part. 
3.1 Evaluation of Vi 

To tame the infrared behavior of the first term in Eq. (2.1), we transform it via the integration- 
by-parts (IBP) relation of Eq. (E.l) into 



Vi 



fob 
p [P 2 f 



n(p)H(p) 



1 



^nn + nn-^-n-/f-n 



d-6% P 2 



1 



1 



with 



n 



- Q [Q 2 ] 2 (P-Q)^ 



(3.1) 



(3.2) 



Using U A = U B + T { p 2 { ^% and U B = f R [H 2 ] 2 {P _ R yi 
are listed in App. B, this is then identically re-written as 



R [R 2 ] 2 (P-R) 2 ■ 

as well as the functions G, I and A that 



5 [Q 2 ] 2 (P - Q)'- 



+ [n - 5 B ] T 



1 



q [Q 2 ] 2 (P - Q) 



+ [n - u A ] 



n-n B 



1 



d-Q 
+ V lb + 



+ UU B + U A U-U A fL B -I^IL-I 2 ^u} (3.3) 
A(^,2,l;0) + \TG(2,l,d)A(^,l,l;2)-0 scaMl J +V la + 



G(l,l,d+1) 

Ad ^ 

2G(l,l,rf+l) -G(2,l,rf+1) 
4d 



A(^,1,1;0)+ G(l,l,d + 1)^4(^,2,1; 2)+ 



+ T G(l, 1, d) A(^, 2, 1; 2)] - Oscaiefrce - i§ A(2, 1, 1; 2) - 7| A(2, 1, 1; 0)} 



7T \ e 1 



(4vr) 4 V4eW T 6e k 



{4^ + 4^ + COnSt + ° (£) }' 



(3.4) 
(3.5) 



- 3 - 



where the const also includes the numerical contributions from V\ a and V±b, which are defined 
in Eqs. (3.6), (3.7) below. For the first term of Eq. (3.3), we have used that 5p tl B = 
2 y^2_ d )/ 2 i s a massless bubble, while for the second term we have used that the 
Q-integration gives a massless bubble, 



To 



[Q 2 ] 2 (P-Qf 



d d q 1 



Sp T 



G(2,l,d) 

rp2\Z-d/2 ' 



leaving exactly solvable 2-loop tadpoles A. The next two (finite) sum-integrals V\ a and V\b 
will be treated numerically in coordinate space, as discussed below. For the fifth and seventh 
terms in Eq. (3.3), we have used that (with U = (1,0)) 



<5p II B = SpgUnUv 



Ru.Ru 



r [r 2 ] 2 (p — Ry 



S Po U^U u {g^A(P 2 ) + P^P U B(P 2 )} = 5 Po A(P 2 



5p 



P Q\iv PfiPy 



dP 2 



RnR U 



5 Po f P 2 R 2 - (PR) 



R [R 2 ] 2 {P - R) 2 dP 2 J R [R 2 ] 2 {P - R) 2 



AdP 2 
$Po 



212 



P 2 (2R 2 - P 2 ) + (P- R) 2 (P + R) 2 - [R 2 ] 
[R 2 ] 2 (P-R) 2 



2R 2 -P 2 _ 8 Po 2G(1, 1, d + 1) - g(2, 1, d + 1) 

Ad J R [R 2 ] 2 (R + P) 2 + scalefree " Id — 3 - d 



(P 2 Y 

while the remaining three terms of Eq. (3.3) are immediately seen to be 2-loop tadpoles A. 
To complete the evaluation of Vi, we need the two (finite) sum-integrals V\ a and V±b, 



1 



T 2 



Vi 6 = ^C^[n-n A ] 



t, [Q 2 ] 2 (P - Q) 2 
n-n B 



(4tt) 4 

T 2 
(4tt) 4 



Via + 0{e) , 
Vib + 0(e). 



(3.6) 
(3.7) 



For Vi a , using the 3d spatial Fourier transforms of [II — 11b] and II' as given in App. A 
(used at Po = 0); integrating over p via Eq. (C.5); letting |r| = x/(2irT) and |r'| = y/(2itT); 
and using Eq. (C.9) for the angular integral: 



V 



la 



oo roo 




1 

36 



dx / dy f(x,0)-f B (x,0) f'(y,0) 
Jo L J 
00 dx r 



x + y — \x — y\ 
xy 



o 



f(x, 0) - f B (x, 0) 3Li 3 (0 - 3C(3) - 2ir 2 x + Qmx 2 + 4x d 



-0.0285014376988(11 



(3.8) 
(3.9) 



where the functions /, fp and /' are listed in Eqs. (A. 10), (A. 11) and (A. 17), respectively. 

For V lb , re-writing 5 Po [II - II A ] = S Po [II - U B ] - ± J x / ^ + 0(e), where 

b = G(l, 1, 3) while the extra integral is unity and introduced here for notational simplicity; 



using the 3d spatial Fourier transforms of [II — n#] and 



n-n B 



(at Pq = 0); integrating 
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over p via Eq. (C.5); letting |r| = x/(2irT) and |r'| = y/(2irT); and using Eq. (C.9) for the 
angular integral: 



o 

00 dx 
x 



fix, 0) - f B (x, 0) - lj [f(y, 0) - f B (y, 0) 
f(x, 0) — /b(x, 0) — 1 x + In (x csch(x)) 



x + y — |x — y| 
xy 



» +1.197038271143294592038702(1) . 
3.2 Evaluation of V2 

This integral is finite, so we write 



(3.10) 

(3.11) 
(3.12) 



v 2 - £i [n - us] [n - n B ] « y 2a + o( e ) . 



(3.13) 



ft [P2] 2 

Using the 3d spatial Fourier transforms of [II— H B ] and [5— IL B ] from App. A; using Eq. (C.4); 
and letting |r| = x/(2irT) and |r'| = y/(2irT), 



V; 



2a 



3 ' 



30 00 2 z-i 



n=l 



-n-y/ x 2 +y 2 +2xyu 



dx T'dyy, 
6 7o % Jo V ' 



1 



f(x,n) - f B (x,n) f(y,n)-f B (y,n) 



71=1 



n(as+i/) n 

1 + nx — ny — e _2ny (l + nx + ny) 



x{ /fori) - / B fon)] f(y,n) - f B (y,n) + xoyj, 



(3.14) 



where in the last step we have integrated over angles via Eq. (C.8), and split the integration 
region into two to avoid taking absolute values. Finally, summing via Eq. (D.l) (which gives 
a large expression containing polylogarithms) and solving the remaining double integral over 
x and y numerically with the help of Mathematica [12], 

V 2a ~ +0.0143560494342(1) . (3.15) 

3.3 Evaluation of V3 

Noting from Eq. (2.3) that Hp is P-independent, we re-write the third term of Eq. (2.1) as 



V 3 



p '^n D [n - ii c ] = u D ^ ij {n - «5 Po n - n c + 5 Po n c } . (3.16) 



Now the first term is a regular 2-loop sum-integral, which reduces via the IBP relation 
Eq. (E.5) to a product of 1-loop tadpoles; the second term is a special 2-loop tadpole A; 
the third term is immediately recognized to be a sum of 1-loop tadpoles by noticing that 
the P-dependence of lie is trivial (powers of Pq and P 2 only); and, for the same reason, the 
fourth term is scale-free and hence vanishes in dimensional regularization. We hence get the 
exact expression 
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+ 



Ad 



((d + l)/| +e - + + 2/3 /?] - O scalcfrec } 



iTe) { 



20 



+ const + O 



to}, 



(An) 4 \ ai T e J 1 e 
where const is somewhat lengthy, so we do not display it here. 



1 /C(3) . C'(-i) 
24 V 5 C(-l) 



7E 



(3.17) 
(3.18) 



3.4 Evaluation of V4 

Expanding [II # — lip] 
well as G(l,l,4-2e) 



■ ( ai T 2 )^} « G(l, 1, 4-2e) e [In + £>( e )] as 



Expanding [11 B - II D ] = G(l, 1, d+l)[{P 2 )' 

^WxWe ^( e2 ))' we see this hitegral is finite, so we write 



v 4 = [n B - u D ] [n - nc] - ^ ^4 + o(e) 



(3.19) 



Using the 3d spatial Fourier transform of [IT — lie*] from App. A; integrating over angles via 
Eq. (C.7); letting |r| = x/(2irT), |p| = \Po\y and [Pol = 2-irTn; and using the exponential- 
integral Eq. (C.3): 



1 00 1 
n 

n=l 



= z In 
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00 
dx 

7-1 



e- 2nx In 



2x a±e 



7-1 



n 167r 2 



Ei(-2nx) } f(x, n) - f c (x, n) 



+ V 4 



la 



16vr 2 
z ln(a 6 
+0.004820184(1) , 



1 00 1 
6 ^ n 

n=l 



00 

d.r 



-2nx 



; ln 



2x 



Ei(-2nx) 



(3.20) 

(3.21) 

f{x,n) - fc(x,n) 

(3.22) 



where the explicit ai-dependence follows from ai-independence of V3 + V4 term with zq 
from Eq. (3.18), and we have introduced an arbitrary parameter oq to parameterize some 
freedom in evaluating the number V 4a . Mathematica [12] actually gives a result for the sum 
for the first term in curly brackets of Eq. (3.20), but none for the Ei term. For obtaining 
the numerical approximation, we have truncated the sum over n at some n max = 30000. In 
order to estimate the magnitude of the remainder, we have interpolated the one-dimensional 
integral in the range n 6 [10 , 10 5 ] with a simple power law n~ a and performed the summation 
from n = 30001 ... 00 analytically, from which we infer that the remainder is of O(10~ 10 ). 

3.5 Evaluation of V5 

In close analogy to the treatment of V3, we identically re-write the fifth term of Eq. (2.1) as 



V5^^^[n-ndn D 



g(i,m+i) 1 _/ (d+i)p$-p 2 



> 



( a2T 2)(3^)/2 M f p pi 



{n-5p u-u c +Sp o u c } . 

(3.23) 
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Now the first term of Eq. (3.23) is a sum of two regular 2-loop sum- integrals, which both 
reduce to zero: first, from the IBP relation of Eq. (E.2) we get 

f p {d+1) p! "^ n = (d+l) 7(211, 20) - /(111, 00) = /(in, 00) , 

while furthermore, the 2-loop sunset sum-integral 7(111,00) vanishes identically via IBP 
identity Eq. (E.3). The second term of Eq. (3.23) is a special 2-loop tadpole A; the third term 
is immediately recognized to be a sum of 1-loop tadpoles by noticing that the P-dependence 
of lie* is trivial (powers of P 2 only); and, for the same reason, the fourth term is scale- free 
and hence vanishes in dimensional regularization. We hence get the exact expression 

V5 = (al^itVi Td{^ ~ °'BP\ - [0 - A ^ 1. 1; 0)1- 

- [G(l, 1, d + 1) ((d + l)ll +t - I° l+e ) + 27? ((d + l)/f - / 2 )] + scalefree } (3.24) 

^(o^) £ {7 + const + ^}' -^(fFiT- ln ^4)' ^ 

where const is somewhat lengthy, so we do not display it here. 
3.6 Evaluation of 



In close analogy to the treatment of V4, expanding [lis— Hd] ~ — — ,\> 2 [In °pT +0(e)] 



we see that the sixth term of Eq. (2.1) is finite, so we write 

Ve = ^'t^ P - Ho] [Sb - S D ] « ^ F 6 + 0(c) . (3.26) 

Using the 3d spatial Fourier transform of [II — Hq] from App. A; integrating over angles via 
Eq. (C.7); letting |r| = x/(2vrT), |p| = \P \y and |P | = 1<KTn: 

roo r oo (o_ rt ,2\ / l(0^\1 s 



1 r°° roo (g yl\ / Ct2 1 (^Ttn) \ g—nx ~ 

V z = ^Y.J dx J dy ry2 + 1) 2 ln [ 2 y 2 + l )yMnxy)—[f{x,n)-fc{x,n) 

(3.27) 

Now, splitting (3 — y 2 ) = 4 — {y 2 + 1) and integrating over y using Eqs. (C.3), (C.2): 

-I 00 poo _ _ 

F 6 = oE / ^n[e- 2nx (L-l)-Ei(-2nx)] f(x,n)-f c (x,n) - 

6 n=l J ° 1 J 

1 00 POO 1 _ , 

" « E / dx - t 6 " 2 ™^ + Ei ("2nx)] /(x, n) - f c (x, n) , (3.28) 
6 ~, Jo x 1 J 



n=l ' 



where L = In (^fgj£). At a 2 = 167r 2 e -T, the last line of Eq. (3.28) equals § » M^£6 ( gee 
Eq. (2.14) of [10], where its numerical evaluation had cost some sweat; see also Eq. (D.27) 
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of [9], where 0.0034814 is quoted). Noting that the Qi2-dependence of Eq. (3.28) follows from 
a^-hrdependence of V5 + Vq, let us re- write it as (with z\ from Eq. (3.25) above) 



V 16V J 

VQa = 3a 3 zx + / dx^e 



-2nx 



3a 3 zi + 



1 



n=l 
00 r 

dx[f(x,0)- f c (x,0) 



36 

f(x,n) - f c {x,n) 

coth(x) — 1 



+ 0.000037813(1) 



(3.29) 



"3 ( . 1 

- n [a 3 + - 

2x V 3 



"3" 



+ a 3 + 



1 



[1 - coth 2 (x)) 



» +0.0025600539026700475010965(1) , 
^66 = -Vea ln(a 4 ) - - / ds;Vn f(x, n) - f c (x, n) 

6J n =l L 

» -0.002619354(1) , 
C = 36(V6a + «i) hi (a 5 ) - 6 / dx - /(x, re) - fc(x, n) 

J ° n=l X L 



(3.30) 



Ei(— 2nx) — e 



-2nx 



hi 



2x \ 
a^n J 



Ei(-2nx) + e 



-2nx 



111 



(3.31 
2x 



+0.0034960718(1) 



(3.32) 



and where we have introduced three arbitrary parameters {03, 04,05} to parameterize some 
freedom we found in evaluating the three numbers {V% a , Veb, C}. Concerning numerical pre- 
cision, see the discussion below Eq. (3.22). As mentioned above, the constant C has already 
been evaluated in [10] (at 05 = 1; see App. B therein). Some further analytic work on these 
constants seems possible (for example, / — fc is n-independent and the corresponding sums 
can be performed; V% a seems to be a suitable candidate for analytic treatment), but the 
numeric result is at this point fully sufficient for our purposes. 

3.7 Evaluation of V7-V9 

The last three terms of Eq. (2.1) are trivial tadpoles, and can immediately be evaluated 
analytically in terms of the functions G and / from App. B as 



V 7 



(4 



^[Uc-U B ]fl B = 2G(l,l,d + l)i? 
p [P z \ 

T 2 ( 1 \ 3e r 1 -1 



d + 1 2 

1 • 



3+e 
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4d 



l 2+e 
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1 T+7E 



96e^ 



V 9 = 



(4tt) 4 V4vrr 2 

, r -Ln s n B = G(i,i,d + i) 2 

p [P z \ 



+ 2 C(-i) 


fC(3) 


96e 




(d+1 2 
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4d ^2+2e 





+ const + OU 



1 l+2e 



(3.33) 

(3.34) 
(3.35) 

(3.36) 
(3.37) 
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M {-48? + + const + 0( £ )}. (3.38) 



(4tt) 4 
4 Result 



Collecting all expressions that have been evaluated in the previous sections in order to re- 
assemble V2 according to Eq. (2.2), setting d = 3 — 2e and expanding, we obtain our final 
result 

^ « (4^4 ^ 96£ 2 [1 + Vie + V22 e 2 + 0(e 3 )] , (4.1) 
V21 = y + 7e + 2^-^, (4.2) 

«^-|V + .(f + af^-sM 2 .) + ^(3,)-|c,3,-|C(3, + 

143 1) <T"(— 1) 

+ -gg- vr 2 + 8 In 2 2 + 16 ln(2vr) + 8 ln(vr) ln(4vr) - 48 7 i + 21 " 6 ((_l) + 



l(Via + V lb ) + V 2a + V 4a + V 6a + Veb + ^C (4.3) 
3 3d 



+ 96 

ps 128.63807196263994701 • • • + 96 [ - 0.370298231(2)] = 93.0894417(2) , (4.4) 



where the Stieltjes constant 71 is defined by £(l+e) = l/e+7E — 7i e+0(e 2 ), and the constant 
V22 contains the seven (sum-)integrals for which we do not have analytic representations. They 
are defined in Eqs. (3.8), (3.11), (3.14), (3.22), (3.30), (3.31) and (3.32), and were evaluated 
numerically in the above. The main uncertainty in the final numerical result Eq. (4.4) is 
dominated by the ones in V 4a and V^. 



5 Conclusions 

We have successfully evaluated a new massless bosonic three-loop sum-integral of mass dimen- 
sion two, which we have named V2, and which represents one of the two remaining unknown 
ingredients for a 3-loop evaluation of the Debye screening mass in hot Yang-Mills theory, 
[1, 8]. 

Our strategy has relied heavily on the particular structure of V2, containing two differ- 
ent one-loop two-point sub- integrals, whose known analytic properties we have repeatedly 
exploited. This particular line of attack has served to evaluate a number of similar (but 
somewhat less involved) cases of sum- integrals in the past. 

We have obtained all divergent terms of V2 analytically, which will be important when 
it is finally used in a concrete physics computation; for the constant part, we had to resort 
to a numerical treatment for some pieces thereof. While it is certainly possible to increase 
the precision of our numerics substantially, we feel that the numbers given above will be 
sufficiently accurate for all practical purposes. It would of course be nice to obtain all terms 
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in V22 analytically, but we have no systematic method to do so for now. In passing, there 
might be some relation of (parts of) V2 to the sum-integral A4±fi that was treated in [9, 10], 
since the constant C of our Eq. (3.32) had already contributed there. However, his might 
also be pure coincidence, and we have not pursued this relation further. 

Finally, we hope that the last remaining sum-integral that is needed at dimension two is 
amenable to similar techniques, and leave its computation for future work. 
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A 3d Fourier transforms 

All (inverse) 3d spatial Fourier transforms derive from Eqs. (22-24) in [7]: 



e -\Q \r ( -L-Y' 1 f s (\Q \ r ) , (A.l) 



W + QIY (STrTiyJ r \\Q 

where f s (m) = ^ — e m K 3/2 „ s (m) , (A.2) 

with special cases of the modified Bessel function of second kind evaluating to fi(m) = 
fzijn) = 1 and fo(m) = f-i{m) = (1 + 1/m). From this one gets a generic formula for e.g. 
(inverse) 3d spatial Fourier transforms of 1-loop IPs: 



Qo _ T(2nT)H2nT^ f d 3 r _ m k 

*-i / - \ i~ l 



J Q 

cm / - \ 1—1 / - \ j — 1 

4(r,\P \) = n fc e-(H+l^|-|Po|)f I \ f . {w) l-J— \ /,-(|P -n|r) . (A.4) 

n=— 00 VII/ V I I / 

Matsubara sums can then be solved using Eq. (19) of [7], as reprinted in Eq. (D.3) below. In 
particular, we get the (inverse) 3d spatial Fourier transforms (see also [10]) 

n<(P) =7^0 I ^e&e-W'fiif, \P \) (A.5) 



(4tt) 2 j 

f(x, n) = coth(x) + n , (A. 6) 

f B (x,n) = n+- , (A.7) 
x 

1 X 

f c (x,n) = n + - + - , (A.8) 
x 3 

n *0P) ~ [ ^e^e-^m \P \) (A.9) 
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f(x,n) = 3coth 3 (x)-3coth(2;)+n(3coth 2 (x)-2)+3n 2 coth(x)+2?i 3 , (A.10) 

/ B (x,n) = 4 + ^ + — + 2n 3 , (A.ll) 

X^ X A X 

f c (x, n ) = — + — + + 2n 3 + n 2 x - - , A.12 

x 6 x l x 5 

f d 3 r 

4(4^ 

f(x,n) = x [coth 2 (x) - 1 + ncoth(x) + n 2 ] , (A. 14) 



fB(x,n) = x 



1 n 9 
-=■ + -+ n -4 



(A.15) 



^on'(P) =5^4 <5p | ^ e ^ e -l^l7~'(r,0) , (A.16) 
f'(x,0) = -xln(l-e- 2x ) . (A.17) 

B Standard integrals 

For convenience, we collect here the functions used above, as defined in [7, 10]. They are the 
zero-temperature massless one-loop propagator 

r( * - ( ^~ d - f ddq 1 - r( ^ ~ Sl)r( ^ ~ g2)r(512 ~ H mu 

<^{si,s 2 ,a) - \p ) J {27r)d[q2]Sl[{q _ p)2]S2 (47r)d /2 r(si)r(s2)r(d _ Sl2) > V- 1 ) 

the one-loop bosonic tadpoles 

J |Q r 2T((2s-a-d) T(s - j ) _^ 1 

S " 7q [Q 2 ] S {2-KTYs-a-d (4vr )rf/2r( S ) ' * - % [Q2]« * ' 1 • j 

and a specific two-loop tadpole 

A(s 1 ,s 2 ,s 3 ) = A(si,s 2 ,s 3 ;0) , (B.3) 



^Oi, s 2 , s 3 ; s 4 ) = X 



pq [Q 2 MP 2 M(^-Q) 2 ] S3 

2^ V°r^ 



2T 2 c(2si 23 - 2d - S4 ) r(si 3 - f )r(«i 2 - f )r(f - si)r( Sl23 - d) 



(B.4) 



(27rT) 2s i23-2d- S 4 (47r) rf r(s 2 )r(s3)r(d/2)r( Sll 23 - d) 

where s abc ... = s c + s b + s c + ... ; 

C Other integrals 

In the main text, we have used the exponential integral, defined as 

je-t (C.l) 
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for the following integrals (the first of which has already been used in Sec. 2.2 of [10]): 

2 ui /"°? V sinfyld) , a . „. , Izla . 

5F 6 /. * VrV 2 lD ^ = - 2| "E.(-2|--|) + 7E + In if . (C.2) 

4 eW ^ J^nfoy, , « = , n No _ . (c . 3) 



tt|z| y (y 2 + l) 2 y 2 + l 2e 1 -T 
Furthermore, we needed one special 3d Fourier transform and its generalization 

d 3 D • 1 e -\ p a\r 

P e ip \ , ■ ' =1-7777, (C.4) 



(2vr)3 (p 2 + J P 2 ) 2 8vr|P | ' 
d 3 D • 1 1 e -l p ol r e -l p ok 2lPnl 



(2vr) 3 (p 2 + P 2 ) n+1 n! ~ p <> 4vrr 4vr (2|P |) 

with / {0 ,i,2,3}(*) = {l/(2x),l,l + x,2 + 2x + 2x 2 /3} , (C.6) 

as well as the following angular averages: 

- [due^ = S -^l, (C.7) 
2 J_ x pr 



- [ due~ n ^ x2+y2+2xyu 



2 J-x 2xyn 



(l+n\x-y\)e~ n ^ x - yl - (l+n(x+y)) e - n( - x+y) ] , (C.8) 



l - fdu , 1 = x+ y-\ x ~y\ . (c.9) 

2 7-i v^ 2 + y 2 + 2x 2/ u 2xy 



D Summation formulae 

When dealing with discrete sums (see also [7]), Zeta functions and polylogarithms enter 
through 

oo oo n 

C (,) = jV', Li ,(x) = ^^ (D.l) 

i i ^ 

n=l n=l 

with Lii(x) = — ln(l — x) (one can use Eq. (D.l) at {s = 1, x = e~ r } and then get the cases 
s < 1 via 3 r ), while hyperbolic functions appear via (m € Z) 

oo 

^ e -(|n|+|n-m|-M)r = | m | + coth ( r ) ) ( D-2 ) 

n=— oo 

00 l m l i p -2|m|r 

£ e-(l"l + l"- m l-l m D>(|n|) = J>(n) + p (_ft r )- ji _ I + e 2|^l-^ ( _ a 2 ,. ) — ^---7-- , (D.3) 

n=— 00 n=0 

where p(x) can be a polynomial. Note that Eq. (D.2) is just a special case of Eq. (D.3). 
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E IBP relations 



We need one 3- loop IBP relation for Sec. 3.1, which can be derived as 
= d pP o Vi 



dVi + 
dV! + 



X P d P] 

JPQR 



5 Po Rl 



pqr p [p 2 ] 2 Q 2 (p-Q) 2 R 2 (p-R) 2 

5 Po R 2 -N 



p QR [P 2 ] 3 Q 2 [(P-Q) 2 ] 2 R 2 [(P-R) 2 ] 2 
5 Po N = 5 Po { - 4pp(P-Q) 2 (P - R) 2 - 2p(p-q)P 2 (P - R) 2 - 2p(p-r)P 2 (P - Q) 2 } 
let pp -»• P 2 , 2p(p-q) -> (P-Q) 2 -Q 2 + P 2 , 2p(p-r) -»• (P-R) 2 -R 2 + P 2 
= 5 Pa P 2 { - 6(P - Q) 2 (P - P) 2 + (Q 2 - P 2 )(P - P) 2 + (P 2 - P 2 )(P - Q) 2 } 
shift Q-^P-Q and R-^P-R 

(«« - 6)V! + £ [Z [ Q2]2( p_ Q)2 n + n|; [ ^ ]2(P _ i?)2 J • 



The 2-loop IBP relations that we need for Sec. 3.3 and Sec. 3.5 can be obtained in the 

[Po] d [Qo] e 
[P 2 ] a lQ 2 ] b [{P~Q) 2 



same way. Letting I(abc,de) = J PQ , p2]a \^}Jup^ 0) i ]c , we can derive e.g. 



= d p po 1(111,00) 

= d/(iii,oo )+ ^^ pP p2Q2( ;_ Q) 



d/(lll,00) + f pQ { [p2]2 ~^_ Q)2 + p^IpIq]^ Mo P ^ Q - P in last term 



d/(lll,00)+^ 

JF 



2 



4p 2 + 2pq 

pq [P 2 \ 2 Q 2 (P-Q) 
r _3("p 2 _ p 2 ) _ pep _ 2a) 

= d/(lll,00) + Y; rp2i2^D ^2 last term odd under P-Q 

jpq [P J y [P - Q) 

= (d - 3) 7(111, 00) +3/(211,20) . (E.2) 

More conveniently, using our computer-algebraic [13] implementation of a Laporta-type [2] 
algorithm to systematically derive and solve such finite-temperature IBP relations [4] , we get 
the special cases 

7(111; 00) = , (for explicit IBP derivation see [6]) (E.3) 

7(211; 20) =0, (E.4) 
1 

d~5 



1(211; 02) = -— — (2/ 1 / 2 + / 2 / 2 °) , (E.5) 



where the second relation confirms Eq. (E.2). 
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